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We present the self-consistent treatment of the simplest, nontrivial, self-gravitating system of degenerate neutrons, protons and 
electrons in /^-equilibrium within relativistic quantum statistics and the Einstein-Maxwell equations. The impossibility of imposing 
the condition of local charge neutrality on such systems is proved, consequently overcoming the traditional Tolman-Oppenheimer- 
Volkoff treatment. We emphasize the crucial role of imposing the constancy of the generalized Fermi energies. A new approach 
based on the coupled system of the general relativistic Thomas-Fermi-Einstein-Maxwell equations is presented and solved. We 
obtain an explicit solution fulfilling global and not local charge neutrality by solving a sophisticated eigenvalue problem of the 
general relativistic Thomas-Fermi equation. The value of the Coulomb potential at the center of the configuration is eV(0) m„c 2 
and the system is intrinsically stable against Coulomb repulsion in the proton component. This approach is necessary, but not 
sufficient, when strong interactions are introduced. 
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1. Introduction 

The insurgence of critical electric fields in the process of 
ravitational collapse leading to vacuum polarization process 
l|] has convinced us of the necessity of critically reexamin- 
ing the gravitational and electrodynamical properties in neutron 
stars. In this light we have recently generalized the Feynman, 
Metropolis and Teller treatment of compressed atoms to the rel- 
ativistic regimes yfl. We have so enforced, self-consistently 
in a relativistic Thomas-Fermi equation, the condition of /3- 



, Ya. B. Zel- 
||3], J. Fer- 
8] for heavy 



equilibrium extending the works of V. S. Popov y [ 
dovich and V. S. Popov [4], A. B. Migdal et al. 
reirinho et al. [7] and R. Ruffini and L. Stella 
nuclei. Thanks to the existence of scaling laws (see [2] and 
|@]) this treatment has been extrapolated to compressed nu- 
clear matter cores of stellar dimensions with mass numbers 
A =! (mpianck/#z n ) 3 ~ 10 57 or M core ~ M Q . Such configurations 
fulfill global but not local charge neutrality. They have elec- 
tric fields on the core surface, increasing for decreasing values 
of the electron Fermi energy Eg reaching values much larger 
than the critical value E c - m\(? j(eK), for Eg - 0. The as- 
sumption of constant distribution of protons at nuclear densi- 
ties simulates, in such a treatment, the confinement due to the 
strong interactions in the case of nuclei and heavy nuclei and 
due to both the gravitational field and the strong interactions in 
the case of nuclear matter cores of stellar sizes. 
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In this article we introduce explicitly the effects of gravita- 
tion by considering a general relativistic system of degenerate 
fermions composed of neutrons, protons and electrons in p 1 - 
equilibrium: this is the simplest nontrivial system in which new 
electrodynamical and general relativistic properties of the equi- 
librium configuration can be clearly and rigorously illustrated. 
We first prove that the condition of local charge neutrality can 
never be implemented since it violates necessary conditions of 
equilibrium at the microphysical scale. We then prove the exis- 
tence of a solution with global, but not local, charge neutral- 
ity by taking into account essential gravito-electrodynamical 
effects. First we recall the constancy of the general relativis- 
tic Fermi energy of each specie pioneered by O. Klein [10]. 
We subsequently introduce the general relativistic Thomas- 
Fermi equations for the three fermion species fulfilling rela- 
tivistic quantum statistics, governed by the Einstein-Maxwell 
equations. The solution of this system of equations presents 
a formidable mathematical challenge in theoretical physics. 
The traditional difficulties encountered in proving the exis- 
tence and unicity of the solution of the Thomas-Fermi equation 

Gil El El El IS El E3l are here enhanced by the neces- 
sity of solving the general relativistic Thomas-Fermi equation 
coupled with the Einstein-Maxwell system of equations. We 
present the general solution for the equilibrium configuration, 
from the center of the star all the way to the border, giving the 
details of the gravitational field, of the electrodynamical field 
as well as of the conserved quantities. 

We illustrate such a solution by selecting a central density 
p(0) = 3.94p nuc , where p nuc 2.7 x 10 14 g crrT 3 is the nu- 
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clear density. We point out the existence near the boundary 
of the core in the equilibrium configuration of three different 
radii, in decreasing order: R e corresponding to the vanishing of 
the Fermi momentum of the electron component; P F = 0, R p 
corresponding to the vanishing of the Fermi momentum of the 
proton component; P F = and R n corresponding to the radius 
at which the Fermi momentum of neutrons vanishes: P F = 0. 
We then give explicit expressions for the proton versus elec- 
tron density ratio and the proton versus neutron density ratio 
for any value of the radial coordinate as well as for the electric 
potential at the center of the configuration. A novel situation 
occurs: the description of the pressure and density is not any- 
longer a local one. Their determination needs prior knowledge 
of the global electrodynamical and gravitational potentials on 
the entire system as well as of the radii R„, R p and R e . This is a 
necessary outcome of the self-consistent solution of the eigen- 
function within general relativistic Thomas-Fermi equation in 
the Einstein-Maxwell background. As expected from the con- 
siderations in yj], the electric potential at the center of the con- 
figuration fulfills eV(0) m n c 2 and the gravitational potential 
1 - e v(0 ^ 2 a; m n lm p . The implementation of the constancy of 
the general relativistic Fermi energy of each particle species 
and the consequent system of equations illustrated here is the 
simplest possible example admitting a rigorous nontrivial solu- 
tion. It will necessarily apply in the case of additional particle 
species and of the inclusion of nuclear interactions: in this cases 
however it is not sufficient and the contribution of nuclear fields 
must be taken into due account. 



2. The impossibility of a solution with local charge neutral- 
ity 

We consider the equilibrium configurations of a degenerate 
gas of neutrons, protons and electrons with total matter energy 
density and total matter pressure 



v=n,p,e 

v 1 2 C' p 2 2 



(1) 



(2) 



where e,(p) = ^c 2 p 2 + mfc 4 is the relativistic single particle 
energy. In addition, we require the condition of /^-equilibrium 
between neutrons, protons and electrons 

jJ„ = H p + fie , (3) 

where P F denotes the Fermi momentum and /j., = 5£/5n, = 
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c 2 (P F ) 2 + m 2 c 4 is the free-chemical potential of particle- 
species with number density n,- = (f^) 3 /(37r 2 fi, 3 ). We now in- 
troduce the extension to general relativity of the Thomas-Fermi 
equilibrium condition on the generalized Fermi energy E F of 
the electron component 



where e is the fundamental charge, V is the Coulomb potential 
of the configuration and we have introduced the metric 



ds 2 = e v(r) c 2 dt 2 - e Mr) dr 2 - r 2 d6 2 - r 2 sin 2 6dip 2 



(5) 



for a spherically symmetric non-rotating neutron star. The 
metric function A is related to the mass M(r) and the electric 
field E(r) = -e^ v+A ^ 2 V' (a prime stands for radial derivative) 
through 



1 



2GM(r) G 



+ - A r 2 E 2 (r) 



(6) 



Thus the equations for the neutron star equilibrium configura- 
tion consist of the following Einstein-Maxwell equations and 
general relativistic Thomas-Fermi equation 
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where a denotes the fine structure constant, V = E F +eV, P em = 
-E 2 /(8tt) and we have used Eq. (|4]i to obtain Eq. dlOt . 

It can be demonstrated that the assumption of the equilibrium 
condition together with the /^-equilibrium condition (O and 
the hydrostatic equilibrium (0 is enough to guarantee the con- 
stancy of the generalized Fermi energy 



E F = z vl2 Hi - mic 2 ■ 



i — n,p,e , 



(11) 



for all particle species separately. Here denotes the particle 
unit charge of the /-species. Indeed, as shown by Olson and 
Bailyn iflil. [l^l . when the fermion nature of the constituents 
and their degeneracy is taken into account, in the configuration 
of minimum energy the generalized Fermi energies E F defined 
by ( fTTT i must be constant over the entire configuration. These 
minimum energy conditions generalize the equilibrium condi- 
tions of Klein [10] and of Kodama and Yamada [20] to the case 
of degenerate multicomponent fluids with particle species with 
non-zero unit charge. 

If one were to assume, as often done in literature, the local 
charge neutrality condition n e (r) = n p (r) instead of assuming 
the equilibrium condition this would lead to V = iden- 
tically (since there will be no electric fields generated by the 
neutral matter distribution) implying via Eqs. (f3) and (0 



E F + E F = e v/2 ( Me +/i p )- (m e + m p )c 2 = E F 
+ (m„ — m e - m p )c 2 — constant. 



(12) 



E e = e v/ n e - m e c - eV — constant , 



(4) 



Thus the neutron Fermi energy would be constant throughout 
the configuration as well as the sum of the proton and elec- 
tron Fermi energies but not the individual Fermi energies of 
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each component. In Fig. Q] we show the results of the Einstein 
equations for a selected value of the central density of a system 
of degenerate neutrons, protons, and electrons in /3-equilibrium 
under the constraint of local charge neutrality. In particular, we 
have plotted the Fermi energy of the particle species in units 
of the pion rest-energy. It can be seen that indeed the Fermi 
energies of the protons and electrons are not constant through- 
out the configuration which would lead to microscopic insta- 
bility. This proves the impossibility of having a self-consistent 
configuration fulfilling the condition of local charge neutrality 
for our system. This result is complementary to the conclu- 
sion of Eq. (4.6) of 0] who found that, at zero temperature, 
only a dust solution with zero particle kinetic energy can sat- 
isfy the condition of local charge neutrality and such a con- 
figuration is clearly unacceptable for an equilibrium state of a 
self-gravitating system. 
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Figure 1: Fermi energies for neutrons, protons and electrons in units of the 
pion rest-energy for a locally neutral configuration with central density p(0) = 
3.94p 

nuc> where p nuc — 2.7 X 10 g cm denotes the nuclear density. 



3. The solution with global charge neutrality 

We turn now to describe the equilibrium configurations ful- 
filling only global charge neutrality. We solve self-consistently 
Eqs. (0 and ([8]) for the metric, Eq. (O for the hydrostatic equi- 
librium of the three degenerate fermions and, in addition, we 
impose Eq. (O for the /^-equilibrium. The crucial equation re- 
lating the proton and the electron distributions is then given 
by the general relativistic Thomas-Fermi equation ( TTOb . The 
boundary conditions are: for Eq. (0 the regularity at the origin: 
M(0) = 0, for Eq. ((9} a given value of the central density, and 
for Eq. ( TTOl l the regularity at the origin n e (0) = n p (0), and a sec- 
ond condition at infinity which results in an eigenvalue problem 
determined by imposing the global charge neutrality conditions 

V{R e )=E F e , y%R e ) = 0, (13) 

at the radius R e of the electron distribution defined by 

Pf(^) = 0, (14) 



2 e v(tf c )/2 _ 2 



= m e c 



2GM(R e ) 2 
1 — m e c . 

c L R e 



(15) 



Then the eigenvalue problem consists in determining the gravi- 
tational potential and the Coulomb potential at the center of the 
configuration that satisfy the conditions (fT3b-([T5b at the bound- 
ary. 

4. Numerical integration of the equilibrium equations 

The solution for the particle densities, the gravitational po- 
tential, the Coulomb potential and the electric field are shown in 
Fig. (f2]i for a configuration with central density p(0) = 3.94p nuc . 
In order to compare our results with those obtained in the case 
of nuclear matter cores of stellar dimensions ||2|] as well as to 
analyze the gravito-electrodynamical stability of the configura- 
tion we have plotted the electric potential in units of the pion 
rest-energy and the gravitational potential in units of the pion- 
to-proton mass ratio. One particular interesting new feature is 
the approach to the boundary of the configuration: three differ- 
ent radii are present corresponding to distinct radii at which the 
individual particle Fermi pressure vanishes. The radius R e for 
the electron component corresponding to P F (R e ) = 0, the ra- 
dius R p for the proton component corresponding to Pp(R p ) = 
and the radius R„ for the neutron component corresponding to 

P*(Rn) = 0. 

The smallest radius R„ is due to the threshold energy for f3- 
decay which occurs at a density ~ 10 7 g crrT 3 . The radius R p 
is larger than R„ because the proton mass is slightly smaller 
than the neutron mass. Instead, R e > R p due to a combined 
effect of the difference between the proton and electron masses 
and the implementation of the global charge neutrality condi- 
tion through the Thomas-Fermi equilibrium conditions. 

For the configuration of Fig. [2] we found R„ 12.735 km, 
R p - 12.863 km and R e R p + 10 3 /l f where A e = h/(m e c) de- 
notes the electron Compton wavelength. We find that the elec- 
tron component follows closely the proton component up to the 
radius R p and neutralizes the configuration at R e without having 
a net charge, contrary to the results e.g in 11911 . 

It can be seen from Fig. [2] that the negative proton gravita- 
tional potential energy is indeed always larger than the positive 
proton electric potential energy. Therefore the configuration is 
stable against Coulomb repulsion. This confirms the results in 
the simplified case analyzed by M. Rotondo et al. in 0], 

From Eq. ( fTTb and the relation between Fermi momentum 
and the particle density P F = (37T 2 h 3 ni) 1 ^ 3 , we obtain the 
proton-to-electron and proton-to-neutron ratio for any value of 
the radial coordinate 
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Figure 2: Top panel: particle number density of neutrons, protons, and electrons 
approaching the boundary of the configuration in units of the nuclear density 
«nuc — 1.6x l0 38 cm~ 3 . Bottom panel: proton and electron Coulomb potentials 
in units of the pion rest-energy eVI(m n c 2 ) and -eV/(m„c 2 ) respectively and the 
proton gravitational potential in units of the pion mass m p <S>lm n where <S> = 

(e'i 2 - 1). 



Figure 3: Top panel: electron number density for r > R p normalized to its 
value at r = R p . Bottom panel: electric field for r > R p normalized to its value 
at r = R p . We have shown also the behavior of the solution of the general 
relativistic Thomas-Fermi equation )10t for two different eigenvalues close to 
the one which gives the globally neutral configuration. 
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where f(r) = (E F + m p c 2 - eV)/(E F + m e c 2 + eV), g(r) = 
(E F + m p c 2 - eV)/(E F + m n c 2 ) and the constant values of the 



generalized Fermi energies are given by 



El = m n c 2 ^ R " )l2 -m n c 2 , 



E F = m^cVW-™ 2 



E F = m e c 2 e v{R ' )/2 - m e c 



m p c + eV(R p ) , 

2 



(18) 
(19) 
(20) 



A novel situation occurs: the determination of the quanti- 
ties given in Eqs. ( [ToT l and ( TT8l necessarily require the prior 
knowledge of the global electrodynamical and gravitational po- 
tential from the center of the configuration all the way out to the 
boundary defined by the radii R e , R p and R„. This necessity is 
an outcome of the solution for the eigenfunction of the general 
relativistic Thomas-Fermi equation (fTOl ). 

From the regularity condition at the center of the star n e (Q) = 
n p (0) together with Eq. ( TToT l we obtain the Coulomb potential 
at the center of the configuration 



eV(0) 



(m p - m e )c 2 



(m p + m e )c 2 
E F + m n c 2 



1 + 



,v(0) 



(m p — m e )c 2 



(21) 



which after some algebraic manipulation and defining the cen- 
tral density in units of the nuclear density rj = p{Q)/p nuc can be 
estimated as 



eV(0) 



m„c 2 ^ 2 



1 + [P F (0)/(m n c)] 2 



(22) 



where we have approximated the gravitational potential at the 



boundary as e 



v(R e )/2 „ e v(R,,)/2 „ p v(R„)/2 



1 . Then for configu- 



rations with central densities larger than the nuclear density we 
necessarily have <?V(0) > 0.35m„-c 2 . In particular, for the con- 
figuration we have exemplified with rj = 3.94 in Fig. |2j from 
the above expression ( l22b we obtain eV(0) - 0.85m ff c 2 . This 
value of the central potential agrees with the one obtained in 
the simplified case of nuclear matter cores with constant proton 
density H. 

5. Conclusions 

We have proved in the first part of this letter that the treat- 
ment generally used for the description of neutron stars adopt- 
ing the condition of local charge neutrality, is not consistent 
with the Einstein-Maxwell equations and microphysical con- 
ditions of equilibrium consistent with quantum statistics (see 
Fig. Q]). We have shown how to construct a self-consistent so- 
lution for a general relativistic system of degenerate neutrons, 
protons and electrons in /^-equilibrium fulfilling global but not 
local charge neutrality. 



Although the mass-radius relation in the simple example con- 
sidered here in our new treatment, differs slightly from the one 
of the traditional approaches, the differences in the electrody- 
namic structure are clearly very large. As is well-known these 
effects can lead to important astrophysical consequences on the 
physics of the gravitational collapse of a neutron star to a black 
hole H. 

Having established in the simplest possible example the new 
set of Einstein-Maxwell and general relativistic Thomas-Fermi 
equations, we now proceed to extend this approach when strong 
interactions are present [21]. The contribution of the strong 
fields to the energy-momentum tensor, to the four-vector cur- 
rent and consequently to the Einstein-Maxwell equations have 
to be taken into account. Clearly in this more general case, the 
conditions introduced in this letter have to be still fulfilled: the 
r-independence of the generalized Fermi energy of electrons 
and the fulfillment of the general relativistic Thomas-Fermi 
equation [21]. In addition, the generalized Fermi energy of pro- 
tons and neutrons will depend on the nuclear interaction fields. 
The fluid of neutrons, protons and electrons in this more gen- 
eral case does not extend all the way to the neutron star surface 
but is confined to the neutron star core endowed with overcriti- 
cal electric fields, in precise analogy with the case of the com- 
pressed nuclear matter core of stellar dimension described in 
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